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We investigate the scaling property of global monopoles in the expanding universe. By directly 
solving the equations of motion for scalar fields, we follow the time development of the number 
density of global monopoles in the radiation dominated (RD) universe and the matter dominated 
(MD) universe. It is confirmed that the global monopole network relaxes into the scaling regime and 
the number per bubble volume is a constant irrespective of the cosmic time. The number density n{t) 
of global monopoles is given by n{t) ~ (0.43±0.07)/t^ during the RD era and n{t) ~ (0.25±0.05)/t^ 
during the MD era. We also examine the peculiar velocity v of global monopoles. For this purpose, 
we establish a method to measure the peculiar velocity by use of only the local quantities of the 
scalar fields. It is found that v ~ (1.0 ± 0.3) during the RD era and v ~ (0.8 ± 0.3) during the MD 

^~N , era. By use of it, a more accurate analytic estimate for the number density of global monopoles is 

f^ ■ obtained. 
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lO \ I- INTRODUCTION 



The grand unified theory based on a simple group predicts magnetic (gauge) monopoles if it breaks to leave the 
U{\) symmetry of the electromagnetism [|l[. Magnetic monopoles are dangerous because they may overdose our 
^^ universe ||l[. However, no magnetic monopoles have been found yet. In fact, the flux of magnetic monopoles is 
fNj severely constrained by cosmological and astrophysical considerations pj. Thus, magnetic monopoles produced in the 

C \ early universe must be diluted away, annihilated, or swept away (the monopole problem). This monopole problem is 
^^ ' one of the motivations of inflation 0] though other solutions are also proposed Ep] . 

^^ \ On the other hand, global monopoles have drawn less attention. However, while magnetic monopoles are dangerous 

^^ for the cosmic history, global monopoles may be favorable because they may produce primordial density fluctuations 

rl^ responsible for the large scale structure formation and the anisotropy of the cosmic microwave background radiation 

0-< (CMB) [||-|l§. Recent observations of the CMB by the Boomerang [|ll| and the MAXIMA ||l2| experiments found the 

O , first acoustic peak with a spherical harmonic multipole I ^ 200 predicted by the standard inflationary scenario. But, 

1^ , they also found a relatively low second peak, which may suggest the contribution of topological defects fl^. Moreover, 

T^ ■ deviations from Gaussianity in CMB are reported in [Q. Thus, though it is improbable for topological defects to 

^ [ become the primary source of primordial density fluctuations, a hybrid model is still attractive, where primordial 

density fluctuations are comprised of adiabatic fluctuations induced by inflation and isocurvature ones induced by 

topological defects. In fact, topological defects can be easily compatible with inflation ||l5Hl7|. 

The key property of global monopoles to contribute primordial density fluctuations properly is scaling, where the 
typical scale of the global monopole network grows in proportion to the horizon scale. [J Then the number density of 
global monopoles is proportional to t^'^ {t : the cosmic time). Here we define the scaling parameter ^ as 

^ = n{t)t\ (1) 

where n(t) is the number density of global monopoles. If ^ becomes a constant irrespective of the cosmic time, we 
can conclude that the global monopole network goes into the scaling regime. 

The mass of a global monopole m is proportional to the distance to the nearest neighborhood antimonopole d 
[m ~ Anri'^d, rj : the absolute magnitude of vacuum expectation values (VEV) of scalar fields]. Roughly speaking, the 
distance is the horizon scale (to be exact d ~ t/S,^/^ ) if the network follows the scaling property. Then the density 
fluctuations produced by global monopoles are given by 



^For both the gauge |lq,h^ and the global string network pot|, the scaling property is confirmed so that density fluctuations 
produced by them become scale invariant. 
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^ (24) T^-^e'^ f ^-) for RD (MD), (2) 
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where ra-p\ = 1.2 x 10^^ GeV is the Plank mass. Thus, the density fluctuations produced by global monopoles become 
scale invariant. One may wonder if the precise value of ^ is not so important because the amplitude of density 
fluctuations depends on only the combination of S^^^^rf. It is true for the case where global monopoles are the 
dominant source of density fluctuations. But that is not the case. As stated earlier, the recent observations indicate 
that the dominant source of density fluctuations is inflation and topological defects contribute to them subdominantly 
In such a case, ^ itself determines the ratio of the above two contributions for a fixed VEV. Thus, the precise value 
of ^ is increasingly important. 

The cosmological evolution of global monopoles was first discussed by Barriola and Vilenkin |q| . They showed that 
the annihilation is so efficient that global monopoles do not overdose our universe unlike magnetic monopoles, but 
that it is not too efficient to survive in our universe. This is mainly because an attractive force works between a 
monopole and an antimonopole. Then, Bennett and Rhie performed the first numerical simulations and found the 
tendency that the number of global monopoles per horizon volume is nearly a constant m . However they used the 
nonlinear a model approximation as equations of motion to evolve the scalar fields. Later, Pen, Spergel, and Turok 
made numerical simulations in both the nonlinear a model approximation and the full potential P]. (See also pO|.) 
However, due to the lack of the computer power, they can run a few realizations so that the scaling property cannot 
be confirmed definitely. Many realizations of numerical simulations are needed to decrease the error and estimate it 
statistically. Furthermore, in order to confirm the scaling property completely, we should pay attention to several 
effects, which may affect the final result, for example, the boundary effect, the grid size effect, the total box size 
dependence, and so on. 

In the previous paper ||21| , we reported the results of our numerical simulations of the global monopole network and 
confirmed that the global monopole network relaxes into the scaling regime both in the RD and the MD universe. In 
this paper, we investigate the cosmological evolution of the global monopole network comprehensively. 

In the next section, we give the formulation and the results of our numerical simulations. From the symmetry 
restoration phase, we follow the evolution of the scalar fields with the 0(3) symmetry, which breaks to generate global 
monopoles. The time development of the number density of global monopoles is examined. Since we need to perform 
a lot of realizations, it is important to establish the method to identify monopoles automatically from the values of 
the scalar fields. We will propose two identification methods and compare the results obtained by both methods. 
We also investigate the peculiar velocity of global monopoles. Our numerical simulations are based on the Eulerian 
view. Therefore, it is very difficult to know where a monopole moves after the long interval enough to measure the 
velocity. We establish the method to measure the velocity of global monopoles by use of only the local quantities of 
scalar fields. In Sec. Ill, we set up the Boltzmann equation for the time development of the number density of global 
monopoles. Using the peculiar velocity obtained from numerical simulations, an analytic estimate for the number 
density of global monopoles is given. In the final section, we give the summary. 

II. NUMERICAL SIMULATIONS 

First of all, we give the formalism of our numerical simulations to follow the evolution of the global monopole 
network. Later we show the results of our numerical simulations and discuss whether the global monopole network 
goes into the scaling regime. Furthermore, the peculiar velocity of global monopoles is investigated. 

We directly solve the equations of motion for scalar fields in the expanding universe, which have the 0(3) symmetry 
at high temperature and later break to generate global monopoles. We consider the following Lagrangian density for 
scalar fields 0°(x)(a = 1,2,3): 

^r] = Ig^.d'^rd''^'' - Kffr,T]. (3) 

Here g^^^, is the fiat Robertson- Walker metric and the effective potential Vcs [4'"' i T] , which represents the typical second 
order phase transition, is given by 



^{4' - vr + ji<y' - v% (4) 



where = ^/(j)'^4>°-, 77 = (J\J\ — {T/TcY and T^ = |vT5(t is the critical temperature. For T > T^, the potential 
T4ff[0°,r] has a minimum at the origin and the 0(3) symmetry is restored. On the other hand, for T < Tc, new 
minima <j> — rj appear and the symmetry is broken, which leads to the formation of global monopoles. 
The equations of motion for the scalar fields 0" in the expanding universe are given by 

cj^-ix) + iHr{x) - -^V20-(x) + A (02(3.) _ ^2) ^a(^) ^ 0^ (5) 

where the dot represents the time derivative and R{t) is the cosmic scale factor. The Hubble parameter H = R{t)/R{t) 
and the cosmic time t are given by 



A5ml^* ' 2H- T^ 



H' = T^^9*T\ t^ — ^^-^ (forRD), 



H^=aiT)^^g.T\ ^^^^^ (for MD), (6) 

with g, to be the total number of degrees of freedom for the relativistic particles. For the MD case, we have 
defined a{T) [a{T) > 1] as a{T) = pmat(r)/piad(r) = ac{Tc/T), where pmat(T) is the contribution to the energy 
density from nonrelativistic particles, Prad(T') is the contribution from relativistic particles at the temperature T, and 
ttc = PraatiTc) / Pi!ui{Tc) ■ Wc also define the dimensionless parameter C as 

^ _ £rd { 45 

<,RD 



a yldni'^g^ 




In our simulation, we take C,rd,md = 10 and 5 to investigate the C dependence on the result. 

We start the numerical simulations from the 0(3) symmetric phase with the temperature T,; = 2Tc, which corre- 
sponds to ti — ic/4(RD) and ti — tc/(2"\/2)(MD). At the initial time {U < tc), we adopt as the initial condition the 
thermal equilibrium state with the mass 



^A(7;2-T2), (8) 

which is the inverse curvature of the potential at the origin at t = ti . 

Hereafter we normalize the scalar field in units of t^ , t and x in units of ti. We set A to be A = 0.25 and normalize 
the scale factor R{t) as i?(l) = 1. Then, the normalized equations of motion for the scalar fields cj)"" are given by 

^a(^) + A^-a(3,) _ iv20-(a;) + A(</)2 - rjlj,)^ = (for RD), 

cp-{x) + -^4>-{x) ~ -Ls7^^-{x) + A(02 _ rjl,^)r - (for MD) (9) 



with mo = ^CflD 7^^471 and 77MI3 = {^/TE/ 12)^ ^mdVT^^TJW^. 

We perform numerical simulations in seven different sets of lattice sizes and lattice spacings in the RD universe and 
the MD universe (See Tables ||and|l|.). In all cases, the time step is taken to be St = 0.01. In the typical case (1), the 
box size is nearly equal to the horizon volume (H^^)^ and the lattice spacing to the typical core size of a monopole 
6x ~ 1.0/(vAcr) at the final time tf. Furthermore, in order to investigate the dependence of C,, we arrange the case 
(7) with C = 5. We have simulated the system from 10 [(2), (3), (5) and (6)] or 50 [(1), (4), and (7)] different thermal 
initial conditions. Also, in order to investigate the effect of the boundary condition (BC), we adopt both the periodic 
BC and the reflective BC [W^(f>°'{x) = on the boundary]. 

A. Number density 

In order to judge whether the global monopole network relaxes into the scaling regime, we follow the time develop- 
ment of £,. If ^ becomes a constant irrespective of the cosmic time, we can conclude that the global monopole network 
goes into the scaling regime. 



First of all, we need to count the number of global monopoles in the simulation box. For the purpose, we must 
establish the identification method of global monopoles because we can obtain only the values of scalar fields. We 
propose two identification methods and later compare the results obtained by using both methods. In the first method 
(I), we use a static spherically-symmetric solution with the topological charge A'^ = 1, which is obtained by solving 
the equation 

dr^ r dr r^ d4> 

with (f>°'{r,6,ip) = (f>{r)x"'/r, x^ = r sin cos (/?, x^ = rsin0sin(^, and x^ — rcosd. The boundary conditions are given 
by 

(/)(r) -^ rj (r ^ oo), 

0(0) = 0. (11) 

One should notice that a point with 0° = for all a's is not necessarily situated at a lattice point. In the worst 
case, a point with <j)°' — lies at the center of a cube. Then, we make the criterion that a lattice is identified with a 
part of a monopole core if the potential energy density there is larger than that corresponding to the field value of 
a static spherically-symmetric solution at r = ^/SSxphys/'^ [(5cCphys = R{t)Sx], that is, the potential energy density at 
the vertices when a static spherically-symmetric monopole lies at the center of cube. Moreover, in order to reduce 
the error, we look on the identified lattices which are connected as one monopole core. In the other method (II), a 
cubic box is regarded as including a monopole if all 0" = (a = 1, 2, 3) surfaces pass through the cubic box, that 
is, the signs of all eight vertices of the box are not identical for each three field (p"". In this method, we also look on 
the identified boxes which are connected as one monopole core. In fact, as shown later, the results with these two 
identification methods coincide very well. 

First of all, we discuss the evolution of global monopoles in the RD universe. The time development of ^rd in the 
cases from (1) to (6) under the periodic BC is described in Figs. ^ Asterisks (*) represent the time development 
of ^uD for the identification method (I). Squares (□) represent the time development of ^hd for the identification 
method (II). As easily seen, the results with two identification methods coincide very well. We also find that after 
some relaxation period, ^hd becomes a constant irrespective of time for all cases. Though all results are consistent 
within the standard deviation, ^no tends to increase as the box size does. This is understood as follows: under the 
periodic BC, a monopole can annihilate with an antimonopole which lies beyond the boundary. Therefore monopoles 
under the periodic BC annihilate more often than those in the real universe, in particular, monopoles annihilate more 
often for smaller box sizes. On the other hand, the time development of ^nu in the cases from (1) to (6) under the 
reflective BC is described in Figs. H. We also find that, (^no tends to become a constant though more relaxation 
period takes. Contrary to the case under the periodic BC, £,iiD tends to decrease as the box size increases. This 
is understood as follows: under the reflective BC, field configurations of a monopole extend with the same phase 
direction beyond the boundary so that a monopole cannot annihilate any antimonopoles which may lie beyond the 
boundary in the real universe. Therefore monopoles under the reflective BC annihilate less often than those in the 
real universe, in particular, monopoles annihilate less often for smaller box sizes. Thus, ^ud takes a larger value in 
a smaller-box simulation due to the boundary effect .n After all, the real number of the monopole per the horizon 
volume lies in between those under the periodic BC and the reflective BC. ^rd of each case is listed in Table ||. From 
the results of the largest-box simulations [case (6)], we can conclude that the global monopole network relaxes into 
scaling regime in the RD universe and ^rd converges to a constant S^rd — (0.43 ± 0.07). We also show the time 
development of S,rd with C = 5 in the case (7) (Fig. |^). S,rd asymptotically becomes a constant (^ro ~ (0.36 ± 0.01) 
under the periodic BC, which is consistent with the above all cases with C = 10 within the standard deviation.^ Hence 
we can also conclude that C, does not change the essential result. 

For the MD case, we also find that after some relaxation period, the number of global monopoles per the horizon 
volume becomes a constant irrespective of the cosmic time under the periodic BC except for the cases (1), (2), and 
(3), in which global monopoles annihilate too much due to the boundary effect. Also, the number of global monopoles 
per the horizon volume becomes a constant irrespective of the cosmic time under the reflective BC except for the 
case (2), where the boundary effect is the most manifest because of the longest time simulation. The tendency of the 



^One may wonder if ^rd increases even after some relaxation period, particular, in the case (2), which is the longest simulation. 
This is also just the boundary effect because the earlier the cosmic time is, the simulation box is larger than the horizon volume. 
^In this case, ^_r_d under the reflective BC also tends to increase due to the boundary effect. 



boundary effect is the same with that for the RD case. Then, from the results of the largest-box simulations [case (6)], 
we conclude that the global monopole network relaxes into scaling regime in the MD universe and S,md converges to 
a constant ^md ~ (0.25 ± 0.05) (see Fig. ^ and ||). We also show the time development of ^md with C^ = 5 in the case 
(7) (Fig. g). £,MD asymptotically becomes a constant £,md — (0.44 ±0.03) under the reflective BC though monopoles 
tend to disappear under the periodic BC due to the boundary effect. This is consistent with the above all cases with 
^ = 10 within the standard deviation. Thus, we have completely confirmed that the global monopole network goes 
into the scaling regime in both the RD universe and MD universe. 

B. Peculiar velocity 

In this subsection, we investigate the peculiar velocity of global monopoles in the expanding universe. In order 
to measure the peculiar velocity, we need to know where a monopole moves at the next step. For the purpose, we 
need to find the method to look on the monopole found at each time as the same. But, generally speaking, it is very 
difficult in case there are a lot of monopoles in the simulation box. 

Then, looking at the matter from another angle, we make best use of the information of scalar fields. Since the 
values and the time derivatives of scalar fields include all informations about monopoles, the velocity of a monopole 
can be represented by only the information of scalar fields, that is, local quantities. In fact, it is possible as shown 
below. First of all, we expand scalar fields (j)°'{x,t) around (j)^{xo,to) up to the first order. 



^"{x, t) ~ 0°(a;o, io) + V(/.°(a;o, to) • {x - xq) + (/)°(a;o, h)[t - to) (a = 1, 2, 3). 



(12) 



The monopole core is identified with the zero of all scalar fields (j>°'. Assuming a monopole lies at xq at the time to, 
the position x of the monopole core at the sufficiently near time t is obtained as the intersection of the following three 
planes, 

A" ■{x~xo) + B''(t-to)^0, (13) 

where A° = W(j>°'{xo, to) and S" = 0''(a;o, to). These equations are easily solved by the Cramer's formula 



{X - Xo)j 

t-to 



Thus, the peculiar velocity of a global monopole v can be estimated as n 
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(15) 



This method has two main sources to generate errors. First of all, our estimate is correct only up to the first order. 
Then, in order to reduce the error due to this approximation, we evaluate the peculiar velocity for the case (3) in all 
the situations because it is the simulation with the highest resolution (the smallest lattice spacing). Next, a monopole 
does not necessarily lie just on the lattice in our simulations. Especially, in the identification method (I), we have 
identified a lattice with a part of the monopole core if the potential energy at the lattice is larger than a critical value. 
Since we have started the simulation from the thermal equilibrium states, there are still small thermal fluctuations 
at late times, which may accidentally lead to the large potential energy for the lattices not corresponding to the 
monopole core. Thus, in the identification method (I), some lattices which have nothing to do with the monopole core 
may be identified with monopole cores. At such lattices, the velocities obtained by the above formula may become 
extraordinarily large, which causes the large error. Then, in order to reduce the errors, we introduce the cutoff for 



This method to measure the velocity of monopoles can apply to that of strings in the same way. In the future publication, 
we will investigate the velocity of the string network with the aid of this method. 



the velocity and abandon the velocities which are larger than a cutoff at estimating the average and dispersion of the 
velocity. We have chosen several cutoff values and investigated their effects on the results. We found that the results 
do not change much if the cutoff is smaller than a value (of course, it need to be larger than unity). Though the 
velocity is at most unity, we set the cutoff to be 1.5 in order to use as many data as possible and reduce the artificial 
effect. 

The time development of the peculiar velocity of global monopoles during the RD era under both the periodic and 
the reflective BCs is depicted in Fig. 0. Though there is still large uncertainty, the peculiar velocity v takes almost 
the same constant asymptotically under both BCs and is given by vrd ~ (1-0 ± 0.3). On the other hand, the time 
development of the peculiar velocity of global monopoles during the MD era is depicted in Fig. H and the peculiar 
velocity v is given by vmd ~ (0.8 ±0.3). Since monopoles disappear at late times under the periodic EC, the peculiar 
velocity is set to zero in such a situation. 

The obtained values of the peculiar velocity are roughly understood as follows. A constant long-range attractive 
force works between a monopole and an antimonopole due to the gradient energy of the scalar fields. Then, the 
monopole is accelerated very much and the relative velocity rapidly gets to the order of unity. In the RD universe, 
the cosmic expansion is not so rapid that the redshift of the velocity due to the cosmic expansion becomes negligible 
and the velocity reaches almost the unity. On the other hand, in the MD era, the universe expands so rapid that the 
velocity is redshifted and it takes a value smaller than the unity. 

III. ANALYTIC ESTIMATE 

In this section, we give a simple analytic estimate for the scaling parameter ^. 

The evolution for the number density of global monopoles n{t) can be described by the following Boltzmann 
equation,!^ 



^ = -P{t)n{t)~3H{t)n{t), 
n(t) 3mn(t) 



(16) 



where R{t) oc i™, P{t) is the probability per unit time that a monopole annihilates with an antimonopole, and T(t) 
is the period it takes for a pair of monopoles at rest with the mean separation l(t) to pair annihilate. The mean 
separation l(i)is given by l{t) = R[t)rs = n{t)~^'^, where r^ is the mean comoving separation. In the previous 



publication |21|, we assumed that the relative velocity between them reaches the order of unity at once because a 
constant attractive force works between a pair of monopoles irrespective of the separation length. In the previous 
section, we have confirmed that the above assumption is basically correct but the peculiar velocity is smaller than 
unity in the MD universe. Then, assuming that the relative velocity is given by the peculiar velocity obtained in the 
previous section and a pair of monopoles does not spiral around each other for a long time, we give a more accurate 
analytic estimate for the number density of global monopoles. The period T{t) is given by the following relation, 

r^+*° dt 



V 



dr, (17) 



/to Rii) ' JO 
where to is the initial time where a pair of monopoles are at rest. Then, the period T{t) reads 

n^)-(^;7M^) (f-*o«r). (18) 

Inserting this result into the Boltzmann equation (|l^), the number density n{t) takes the following asymptotic value: 

33(l-m) (i_^)3(2-m) 



^A similar discussion was done in ftl,E2|. 



From the above asymptotic form, first of all, we find that the number density n{t) is proportional to the inverse of the 
cosmic time cubed, t~^ , which implies that ^ becomes a constant irrespective of the cosmic time. ^ is also estimated 
as 

3 3(l-m) Q _ \3{2-m) 

6 U^n^ 

Inserting m ~ 1/2 and vrd ^ (1.0 ± 0.3), £,rd ^ (0.45 ± 0.22). On the other hand, inserting m — 2/3 and 
vmd ^ (0.8 ±0.3), £,MD ~ (0.17 ±0.13). Thus, (, obtained by the analytic estimates can well reproduce that obtained 
from the numerical simulations both in the RD universe and the MD universe. 

IV. SUMMARY 

In this paper, we have discussed the evolution of the global monopole network in the expanding universe. We have 
completely confirmed that the global monopole network relaxes into the scaling regime, where the number of global 
monopoles per the horizon volume is a constant. The scaling parameter ^ is given by ^rd — (0.43 ± 0.07) in the RD 
universe and £,md — (0.25 ±0.05) in the MD universe. We also investigated the peculiar velocity of global monopoles. 
First of all, we established the method to measure the peculiar velocity by using only the local quantities of the scalar 
fields. This method compensates the weak point of the Eulerian view which our numerical simulations are based on, 
that is, we cannot follow the motion of each monopole in detail. We find that the peculiar velocity also becomes a 
constant irrespective of the cosmic time and is given by vrd ^ (1-0 ± 0.3) and vmd ^ (0.8 ± 0.3) though there is still 
large uncertainty. By use of the Boltzmann equation for the time development of the number density and the peculiar 
velocity obtained from numerical simulations, we give a simple analytic estimate for the number density, which can 
well reproduce the results from the numerical simulations up to the proportional coefficient ^. 
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FIG. 1. The time development of ^hd in the cases from (1) to (6) under the periodic BC for the RD case. Asterisks (*) 
represent time development of ^no for the identification method (I). Squares (□) represent time development of ^jiu for the 
identification method (II). The vertical lines denote a standard deviation over difi'erent initial conditions for the identification 
method (II). 
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FIG. 3. The time development of S^rd in the RD universe for the case (7) under the periodic BC and the reflective BC. 
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FIG. 7. The time development of the pecuhar velocity in the RD universe for the case (3) under the periodic and the reflective 
BC. Asterisks (*) represent time development of the peculiar velocity v for the identification method (I). Squares (□) represent 
time development of the peculiar velocity v for the identification method (II). The vertical lines denote a standard deviation 
over different initial conditions for the identification method (II). 
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FIG. 8. The time development of the pecuhar velocity in the MD universe for the case (3) under the periodic and the 
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TABLE I. Seven different sets of the simulations for the RD case. 



Case 


Lattice 


Lattice spacing (Sx) 


c 


Reahzation 


Box svLc/ H ^ 


C 


i 




number 


[unit = tiR{t)] 






(at final time) 


(periodic B.C.) 


(reflective B.C.) 


(1) 


128^ 


^/3/10 


10 


50 


l(at 125) 


0.28 ±0.19 


1.00 ± 035 


(2) 


256^ 


^/6/20 


10 


10 


l(at 250) 


0.31 ±0.20 


0.83 ±0.15 


(3) 


256^ 


\/3/20 


10 


10 


l(at 125) 


0.35 ±0.21 


0.71 ±0.33 


(4) 


128^ 


^/5 


10 


50 


2(at 125) 


0.37 ±0.06 


0.71 ±0.11 


(5) 


256^ 


V3/10 


10 


10 


2(at 125) 


0.36 ± 0.07 


0.58 ±0.12 


(6) 


256^ 


^/3/5 


10 


10 


4(at 125) 


0.36 ±0.01 


0.50 ±0.03 


(7) 


128'^ 


V6/W 


5 


50 


l(at 250) 


0.36 ±0.17 


0.61 ±0.21 



TABLE II. Seven different sets of the simulations for the MD case. 



Case 


Lattice 


Lattice spacing {5x) 


C 


Realization 


Box size/H ^ 


C 


c 




number 


[unit = tiR{t)] 






(at final time) 


(periodic B.C.) 


(reflective B.C.) 


(1) 


128=^ 


3(100)'/V256 


10 


50 


l(at 100) 


Disappearance 


0.75 ±0.36 


(2) 


256^ 


3(200)'/Vl28 


10 


10 


l(at 200) 


Disappearance 


0.77 ±0.10 


(3) 


256^ 


3(100)^/V512 


10 


10 


l(at 100) 


Disappearance 


0.82 ±0.50 


(4) 


128^ 


3(100)^/Vl28 


10 


50 


2(at 100) 


0.19 ±0.09 


0.42 ±0.10 


(5) 


256^ 


3(100)^/^256 


10 


10 


2(at 100) 


0.21 ±0.09 


0.41 ± 0.05 


(6) 


256^ 


3(100)^/^128 


10 


10 


4(at 100) 


0.20 ±0.02 


0.30 ±0.02 


(7) 


128=^ 


3(200)^/^256 


5 


50 


l(at 200) 


Disappearance 


0.44 ± 0.03 
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